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From open quantum walks to unitary quantum walks 
Chaobin Liu Q 
Abstract 

We present an idea to convert to a unitary quantum walk any open 
quantum walk which is defined on lattices as well as on finite graphs. 

This approach generalizes to the domain of open quantum walks (or 
quantum Markov chains) the framework introduced by Szegedy for 
quantizing Markov chains. For the unitary quantum walks formulated 
in this article, we define the probability and the mean probability of 
finding the walk at a node, then derive the asymptotic mean proba¬ 
bility. 

keywords: Quantum Markov chains, open quantum walks, quantum 
walks. 


1 Introduction 

Markov chains or random walks (Markov chains on lattices or graphs), as 
statistical models of real-world processes, have broad applications in various 
holds of mathematics, computer science, physics and the natural sciences. 
An exact quantum extension of Markov chains (named quantum Markov 
chains ) was recently presented by Gudder [T] , at which the author defines 
a transition operation matrix as a matrix whose entries are completely pos¬ 
itive maps whose column sums form a quantum operation. More recently, 
Attal et al. PIE] introduced a formalism for discrete time open quantum 
walks (OQW), which are formulated as quantum Markov chains on lattices 
or graphs. For a short introduction to OQW and some of the recent devel¬ 
opments on OQW (quantum Markov chains), readers may refer to [3]. 

It is noteworthy that the concept of quantum walks, as a quantum coun¬ 
terpart of random walks, has been proposed and well developed prior to 
the introduction of OQW. For a lively and informative elaboration of the 
history of such quantum walks and their connection to modern sciences, the 
reader is referred to [61 !Zl El El dnum H21HH HH H5l ng and references cited 
therein. 

Random walks (or Markov chains), in particular, can be used in computer 
science to formulate search algorithm. To invent a more efficient search 
algorithm, Szegedy [5] developed a generic method for quantizing classical 
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algorithms based on random walks. He shows that under certain conditions, 
the quantum version gives rise to a quadratic speed-up. Enlightened by 
the techniques (namely, quantizing a random walk to generate a unitary 
quantum walk) employed by Szegedy, it is possible, in the present context, 
to extend it to the domain of OQW. We will show how to convert open 
quantum walks to unitary quantum walks. 

The structure of this article is outlined as follow: In section 2, we briefly 
review the formalism of OQW, then we show how to convert an open quan¬ 
tum walk to a unitary quantum walk in section 3. Under the framework 
provided in section 3, we define the probability and the mean probability of 
the walk at nodes and discuss asymptotic distributions of the unitary quan¬ 
tum walks in section 4. In the following section, we show how to recover 
Szegedy’s quantum walks by our approach. Finally we summarize and offer 
some remarks on the notions of Markov chains, open quantum walks and 
unitary quantum walks. 


2 Review of open quantum walks 

Before we proceed to show how to quantize an open quantum walk, let 
us to review briefly the formalism of open quantum walks introduced by 
Attal et al. EE]. One considers OQW on a directed graph G(V,E). V 
is the set of vertices of the graph G, E is the set of oriented edge of G, 
E = {( j,k ) : j,k G V}. The position space of states corresponding to 
the dynamics on the graph is denoted by Hy = C^l with the standard 
orthonormal basis {| j)}jev- The set of all linear operators on the Hilbert 
space Hy is denoted by 23 ("Hy)- The set of degrees of freedom (or called 
“coin”) is denoted by C = {ci, C 2 ,..., c n }, the set of all linear operators on 
the Hilbert space He = C n = span{|ci), |C 2 ), |c n )} is denoted by 23("Hc)- 

Let T>(Hc) C *B(Hc) denote the set of positive operators p : C n —» C n 
with Tr(p) = 1. The operators p are the so-called density operators. Thus, an 
overall state of the quantum walker can be described on the tensor product 
of the Hilbert spaces: 23 (He) <8> IB (Hv )• 

To describe the dynamics of the quantum walker, for each edge (j,k), 
one introduces a bounded operator Bj G 23 (He)- This operator describe 
the change in the internal degree of freedom of the walker due to the shift 
from vertex j to vertex k. By imposing for each vertex j that 

V Bf B* = (1) 

k 

One makes sure, that for each vertex of the graph j G V there is a corre¬ 
sponding completely positive map (quantum operation) on the positive oper¬ 
ators of 23 {He) ■ M.j{r) = Yhk ■ Since the operators B* act only on 

*B(Hc) but don’t perform transitions from vertex j to vertex k, and operator 
M* G 23 (He) ® 23 (Hv) is introduced with the form M* = Bj <g) \k)(j\. It 


is evident that, with the Eq. © being satisfied, one has yE k = I. 

This condition defines a completely positive map for a density operator 
pe<B(Hc)®&(H v ), he., 

Mj pMj. 

k j 

The above map defines the transition matrix for the open quantum walk. 
Given p 0 £ ( B(Hc)®%$(Hv), where tr(p 0 ) = 1, the expression p t = Ai t po 
is called the state of the walker at time t. The corresponding OQW with 
initial state p 0 is represented by the sequence {pt}£fi 0 . 

3 How to convert an open quantum walk to 
a unitary quantum walk 

In what follows, we proceed to “quantize” an open quantum walk. For any 
completely positive map M. defined above, we will define a corresponding 
unitary quantum walk, residing on the augmented Hilbert space 23 (He) ® 
Hv <8> Hv with the Hilbert-Schmidt inner product, defined by 

(A, B) =tr (A'B) 


To define this walk, we first introduce the states: 

IV’j) B f B j ® 1 3) ® I*), (2) 

* k 

For j = l,2,...,|Vj. Each such state is normalized, to see this, we calcu¬ 
late the norm of |^-). |||^)|| = \/(^#j» = yj B f B j) = 1 by the 

condition given in Eq. (1). Then we define 

M 

(3) 

j 

which is the projection on span{|^-) : j — 1, 2,..., |Vj}, denoted by "Hp. This 
is a subspace of the augmented Hilbert space 23 (He) <8> Hv <S> Hv- Finally 
let us define 


\v\ 

S :=I n ® ^ \j,k)(k,j\ (4) 

j,k =1 

as the operator that swaps the two registers. 

With the operators defined above, a single step of the quantum walk is 
defined as the unitary operator U := 5(2n — 1). Given |a 0 ) G *B(Hc) <S> 





'Hy®'Hy, where |||a 0 )|| = 1, the expression \a t ) = U^olq) is called the state 
for the walk at time t. The corresponding quantum walk with initial state 
|a 0 ) is represented by the sequence {|cq)}£h 0 . 

In order to understand the behavior of a quantum walk, one needs to 
know the spectral properties of the unitary operator U. It will be helpful 
to begin with the study of an \V\ x |D| matrix D = (djk), as a linear 
transformation on the space T~Ly = C^ v L The entries of this matrix is defined 
as follows: 


= < 5 > 

Let us then define an operator A from the space 'Hy to Tty: 

|V"( 

A = ^2\^j){J\ ( 6 ) 

j = 1 

The following identities describe the relationships among these operators: 

A f A = I, AA f = n, A^SA = D 

Since D is symmetric by its definition, without loss of the generality, we 
may assume that, via the Spectral Decomposition, D = Y r A r |uy)(uy| + 
Y s K)( M s| - Yt I v t)( v t) where A r e (-1,1), and (|uy), \u s ), |u t )} is an or¬ 
thonormal basis for Hy. 

Due to that = S\i/)j) and US\ij}j) = 2 Yk djkS\ipk) — \ipj), we see 

that the subspace 'H i ys = span{|^-), S\if)j) : j G V} is invariant under U. 
Notice that Y r A\w r )(w r \Al + Y S A\u s )(u s \Al + Yt M v t){ v t\^ = AA f = II, 
the subspaces span{A| w r ), A|w s ), A\v t )} = span{|^-) : j G V} and thus 
= span{A|uy), SA\w r ),A\u s ), SA\u s ), A\v t ), SA\v t )} 

It can be shown with not much difficulty that the following assertions 
hold: 

1. A\w)-e ±l ™ ccosX SA\w) is an eigenvector of U with corresponding eigen¬ 
value e ±iarccosA . 

2. A\u) = SA\u), and A\u) is an eigenvector of U with corresponding 
eigenvalue 1. 

3. A\v) = — S34|u), and A\v) is an eigenvector of U with corresponding 
eigenvalue —1. 

The last two facts imply that — span{A|uy), SA\w r ), A\u s ), A\v t )}. 

Since { A\w r ) — e ±,iaiccosXr SA\w r }, A\u s ), A\v t )} forms an orthogonal set, 
= span{A|uy) — e ±iajTCCOsXr SA\w r ),A\u s ),A\v t )}. After normalizing 



each vector in this orthogonal set, we obtain an orthonormal basis for the in¬ 
variant subspace which is given by the set {A\w+), A\w~), A|w s ), A\v t )}. 
Here 

1. A\w+) = (. A\w r ) - e ia ’ IccosXr SA\w r )/ y /2-2X 2 r 

2. A\w~) = ( A\w r ) - e —Grecos A r SA\w r ))/y/2 ~ 2A ? 2 

To summarize the arguments made above, H^^s — spa.n{), S\i/jj) : j G V} 
as an invariant subspace under U, can be recasted by 

span{A|w+), A\w~), A\u a ), A\v t }} 

where this spanning set is the collection of the orthonormal eigenvectors of 
U associated with the key operator D. 

Let us decompose the Hilbert space 93 (' He ) <8> Hv ® Hv into H^.s and 
its orthogonal complement H^; s , i.e., ^8 (He) <E> Hv ® Hv = H^.s © H^ s . 
It is easy to check that the actions of U and U 2 on H^ s are —S (thus 
H^ s is invariant under U) and the identity I, respectively. Therefore, the 
nontrivial action of U only takes place on the subspace H^.s of a dimension 
less than or equal to 2\V\ (this maximum dimension can be achieved only if 
D does not have both 1 and —1 as its eigenvalues). Based on the aforesaid 
observation, we may confine the initial state of the quantum walk to the 
subspace H^,Si which is spanned by the set of the orthonormal eigenvectors 
of U: {A\w+’),A\w~),A\u s ),A\v t )}. 

4 Asymptotic distribution of the quantum 
walks 

We now turn to study the evolution of the quantum walk. Beginning with 
the initial state |ao), the state of the unitary quantum walk at time t is 
| a t ) = U^ao). Since U is unitary, in general the limit linp^oo \a t ) does not 
exist. Now consider instead the probability distribution on the states of the 
underlying open quantum walks induced by | a t ), 

Definition. P t (j\a 0 ) = y) fc (I n 2 <g> \j,k)(j,k\a t )(a t \). Here P t (j\a 0 ) is the 
probability of finding the walk at the node Vj at time t. 

As a matter of fact, P t usually do not converge either. However, the 
average of P t over time is convergent. We define: 

Definition. Pr(j|«o) = y Y^t=i Pt(j\ a o)- This is the mean probability of 
finding the walk at the node v 3 over time interval [1, T], 

For the sake of brevity, we denote the set of eigenvalues of U by {4>i}, 
and the set of the corresponding eigenvalues of U by {/p}- Via a routine 
reasoning, we can arrive at a theorem regarding the asymptotic distribution 
of the unitary quantum walks. 



Theorem Given an open quantum walk on the state space V with the 
transition matrix A4, the induced quantum walk is defined as | a t ) = U^olq) 
where the initial state |a 0 ) — '52i(4 > i\° l o) |0z)> then 

lim P T {j\ao) = y'y'(^|«o)(Qio|0m)(In® \j,k)(j,k\,\<j>i)(<j> m \) 

T\-too zJ z ' 

k l,m 

where the first sum is over all values of k, and the second sum is only on 
pairs l,m such that pi = p m . 


Example Let us consider a simple open quantum walk on the graph with 
two vertices (see Figured]), its transition operator is given by 
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( 8 ) 
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( 10 ) 

( 11 ) 


( 12 ) 


Here k — 1,2,.... Therefore p^ does not exist in this case! However, if one 
performs measurements of the position of the walker at each node, the quan¬ 
tum trajectories [|3] of the walker will converge to Poo( v i) — CPoo^) = 0 in 
each case. These quantum trajectories imitate the corresponding classical 
Markov chain (illustrated in Figure [2]) with transition matrix given by 


P 


1 1 
0 0 


(13) 


Under the framework described in section 3, we will convert this OQW 
to a unitary quantum walk (see Figure ED- 
By Eq .©, 
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Figure 1: Open quantum walk 



Figure 2: Markov chain 



Figure 3: Quantum walk 
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|2,1), (14) 


Case 1. We choose |a 0 ) = ^(|Vh) + IV^)) = \\ ® 11,1) + §I 2 ® |2,1), 


then we have 


|«i) = U\a 0 ) = |l 2 ® |1,1) + fl 2 ® |1, 2 ), and Pi(l|a 0 ) = 1, Pi( 2 |a 0 ) = 0. 
\a 2 ) = U\ai) = \l 2 ® |1,1) - \l 2 ® 1 2 , 1 ), and P 2 (l|a 0 ) = P 2 (2|a 0 ) = \- 
ja 3 ) = U\a 2 ) = ll 2 ® |1,1) - 5 I 2 ® |1,2), and P 3 (l|a 0 ) = l,P 3 (2|a 0 ) = 0. 
|a 4 ) = U\a 3 ) = |I 2 ® |1,1) + \I 2 < 8 > | 2 , 1 ), and P 4 (l|a 0 ) = P 4 (2|a 0 ) = 


It is seen that the states of the quantum walker {cq}£h 0 are periodic, and 

the period is 4. The limiting distribution is: Poo(l|ao) = f and P oo (2|a 0 ) = 
1 

4' 

Case 2. We choose |a 0 ) = ^H -2 ® 11,2), then we have 

|«i) = U\a 0 ) = ® |2,1), and Pi(l|a 0 ) = 0, Pi(2|« 0 ) = 1. 

\a 2 ) = U\ai) = —^I 2 ® 11, 2), and P 2 (l|a 0 ) = 1, P 2 (2|ai) = 0. 

|a 3 ) = U\a 2 ) = ® |2,1), and P 3 (l|a 0 ) = 0, P 3 (2|a 0 ) = 1. 

|a 4 ) = U\a 3 ) = -^I 2 ® 11, 2), and P 4 (l|a 0 ) = 1, Pi(2|ao) = 0. 


Again the the states of the quantum walker are periodic, and the period 
is also equal to 4. However, the limiting distribution is different from the 
one in case 1 . In this case, it is: P oo (l|a 0 ) = -Poo(2|a 0 ) = \- 

In summary, the mean probability distributions of the quantum walks 
are convergent by Theorem 1. The asymptotic mean distributions are de¬ 
pendent on initial states because the underlying network is not connected in 
the example we discussed. In general, it is speculated that the asymptotic 
distribution for the quantum walk may be unique if the underlying network 
is strongly connected. 


5 Recovering Szegedy’s quantum walks 

Let us recall how Szegedy [5j quantizes a random walk (Markov chain) to 
produce a unitary quantum walk (we call it Szegedy’s quantum walk). 

A discrete-time random walk on an A r -vertex graph can be represented 
by an N x N matrix P in which the entry p^j represents the probability 
of making a transition to k from j. Let u be the probability vector which 
represents the starting distribution. Then the probability distribution after 
one step of the walk becomes Pu. To preserve normalization, we must have 
J2k=i Pkj — 1 ) suc h a matrix is often said to be stochastic. 

For the N x N stochastic matrix P, one can define a corresponding 
discrete-time quantum walk by a unitary operation on the Hilbert space 
C N < 8 ) C N . To define this unitary operation, one first introduces the states 






N 

Hj) : = Y1 VPkj\j> k ) 

k =1 

for j = 1,2,...,N. Each such state is normalized because P is stochastic. 
Then one defines 


N 

3 = 1 

which is the projection onto span{^) : j — 1, 2, N}, and finally one lets 

N 

s ■■= \p k )( k >j\ 

j,k=l 

be the operator that swaps the two registers. Then a single step of Szegedy’s 
quantum walk is defined as the unitary operator U : = S'(2fl — 1). 

To see how the operator U acts on the underlying Hilbert space Py <8> 
Py = span{|j, k) : j, k G V}, it suffices to check U(\jo,ko)). By the dehni- 
tion for U defined above, we have 

U(\j 0 ,k Q )) = 2 y/p kojo VP*h\ k iio) ~ \ko,jo) ( 15 ) 

k 

Now we turn to the setup of converting an open quantum walks to a 
unitary quantum walks presented in this article. 

Consider the stochastic transition matrix P = (pkj), and any family of 
unitary operators Uj on C n for j, k G V, we set 



ft is seen that Y^k=\ B'f B* = l n for all j. 

Employing the unitary operation U defined in section 3, it can be calcu¬ 
lated that 

U(I n ® \jo, k 0 }) = In® [2y/Pkoj 0 ^2 y/Pkfr \k,jo) - 1^0, Jo)] (16) 

k 

From Eqs. f|T5]l and (TlGl) . we can see that the two unitary operators are 
identical over %v ® Pv- It is noted that the norrned spaces Py <S> Pv (the 
underlying space for Szegedy’s walk) and span{I n } <g) Py <g> Py (a subspace 
of the augmented Hilbert space 23 (Pc) Py <8> Py) are isometrically iso¬ 
morphic, therefore Szegedy’s quantum walk can be considered as a special 
case of the quantum walk defined in this article. 




6 Summery and some remarks 

We present an scheme to convert an open quantum walk to a unitary quan¬ 
tum walk. This approach extends to the domain of open quantum walks 
(or quantum Markov chains) the framework introduced by Szegedy [5] for 
quantizing Markov chains. Open quantum walks can be viewed as an exact 
quantum extension of random walks (Markov chains). To illustrate the re¬ 
lationships among the notions of random walks, quantum walks and open 
quantum walks, we draw a diagram below (the work done in this article is 
marked in grey). 



For the quantum walks we introduced in this article, we define the prob¬ 
ability and the mean probability of finding the walker at a node, then we 
obtain a theorem regarding the asymptotic mean probability distribution. 


References 

[1] S. Gudder, Quantum Markov chains , J Math. Phys., 49, 072105(2008). 

[2] S. Attal, F. Petruccione and I. Sinayskiy(2012), Open quantum walks on 
graphs, Phys. Lett. A 376 1545. 






[3] S. Attal, F. Petruccione, C. Sabot and I. Sinayskiy(2012), Open quantum 
random walks, J. Stat. Phys. 147 832. 

[4] I. Sinayskiy and F. Petruccione (2013), Open Quantum Walks: a short 
introduction J. Phys.: Conf. Ser. 442 012003. 

Y. Feng, N. Yu and M. Ying, checking quantum Markov chains, Journal 
of Computer and System Sciences Volume 79, Issue 7, November 2013, 
Pages 11811198. 

C. F. Lardizabal and R. R. Souza, On a class of quantum channels, open 
random walks and recurrence , arXiv:1402.0483. 

R. Carbone and Y. Pautrat, Homogeneous Open Quantum Random 
Walks on a lattice, arXiv:1408.1113. 

R. Carbone and Y. Pautrat, Open Quantum Random Walks: reducibility, 
period , ergodic properties , arXiv:1405.2214, 

R. Sweke, I. Sinayskiy and F. Petruccione, Dissipative preparation of 
large W states in optical cavities, Phys. Rev. A 87, 042323. 

N. Konno and H. J. Yoo, Limit Theorems for Open Quantum Random 
Walks, Journal of Statistical Physics January 2013, Volume 150, Issue 2, 
pp 299-319. 

I. Sinayskiy and F. Petruccione, Properties of open quantum walks on Z, 
Phys. Scr. 2012 014077. 

I. Sinayskiy and F. Petruccione, Microscopic Derivation of Open Quan¬ 
tum Walk on Two-Node Graph, Open Syst. Inf. Dyn. 20, 1340007 (2013). 

I. Sinayskiy and F. Petruccione, Efficiency of open quantum walk imple¬ 
mentation of dissipative quantum computing algorithms, Quantum Infor¬ 
mation Processing October 2012, Volume 11, Issue 5, pp 1301-1309. 

M. Bauer, D. Bernard and A. Tilloy, Open quantum random walks: Bista¬ 
bility on pure states and ballistically induced diffusion, Phys. Rev. A 88, 
062340 (2013). 

M. Bauer, D. Bernard and A. Tilloy, The Open Quantum Brownian Mo¬ 
tion , arXiv:1312.1600. 

C. Pellegrini, Continuous Time Open Quantum Random Walks and Non- 
Markovian Lindblad Master Equations, Journal of Statistical Physics, 
February 2014, Volume 154, Issue 3, pp 838-865. 

R. Sweke, I. Sinayskiy and F. Petruccione, Dissipative preparation of 
generalized Bell states, 2013 J. Phys. B: At. Mol. Opt. Phys. 46 104004. 

C. Ampadu, Return Probability of the Open Quantum Random Walk with 
Time-Dependence, Connnun. Tlieor. Phys. 59 (2013) 563567. 

P. Sadowski and L. Pawela, Central limit theorem for reducible and irre¬ 
ducible open quantum walks, arXiv:1412.5086, 


C. Ampadu, Averaging in SU(2) open quantum random walk , 2014 Chi¬ 
nese Phys. B 23 030302. 

S. Xiong and W.-S. Yang, Open Quantum Random Walks with Decoher¬ 
ence on Coins with n Degrees of Freedom, Journal of Statistical Physics 
August 2013, Volume 152, Issue 3, pp 473-492. 

I. Sinayskiy and F. Petruccione, Quantum optical implementation of open 
quantum walks, Int. J. Quantum Inform. 12, 1461010 (2014). 

M. Schuld, I. Sinayskiy and F. Petruccione, Quantum walks on graphs 
representing the firing patterns of a quantum neural network, Phys. Rev. 
A 89, 032333 (2014). 

S. Attal, N. Guillotin-Plantard and C. Sabot, Central Limit Theorems 
for Open Quantum Random Walks and Quantum Measurement Records, 
Annales Henri Poincare, January 2015, Volume 16, Issue 1, pp 15-43. 

[5] M. Szegedy, Quantum Speed-Up of Markov Chain Based Algorithms, 
Foundations of Computer Science, 2004. Proceedings of the 45th Annual 
IEEE Symposium on Foundations of Computer Science, Pages 32-41. 

[6] A. Ambainis, E. Bach, A. Nayak, A. Vishwanath and J. Watrous, One¬ 
dimensional quantum walks, in Proceedings of the thirty-third annual 
ACM symposium on Theory of computing, STOC ’01 (ACM, New York, 
NY, USA, 2001) pp. 37-49. 

[7] J. Watrous, Quantum simulations of classical random walks and undi¬ 
rected graph connectivity, Journal of Computer and System Sciences 62, 
pp. 376391 (2001), eprint cs.CC/9812012. 

[8] D. Aharanov, A. Ambainis, J. Kempe and U. Vazirani, Quantum walks 
on graphs, in Proceedings of the 33rd Annual ACM Symposium on The¬ 
ory of Computing, (ACM, New York, 2001), pp.50-59. 

[9] E. Far hi and S. Gutmann, Quantum computation and decision trees, 
Phys. Rev. A 58, 915928 (1998). 

[10] J. Kempe (2003), Quantum random walks - an introductory overview, 
Contemp. Phys. 44, 307. 

[11] A. M. Childs, R. Cleve, E. Deotto, E. Farhi, S. Gutmann and D. A. 
Spielman, 2003, Exponential algorithmic speedup by quantum walk, Pro¬ 
ceedings of the 35th ACM Symposium on Theory of Computing, pp. 
5968, eprint quant-ph/0209131. 

[12] N. Konno, A new type of limit theorems for the one-dimensional quan¬ 
tum random walk, Journal of the Mathematical Society of Japan, 57: 
1179-1195 (2005). 

[13] V. Kendon, Decoherence in quantum walks - a review, Struct, in Comp. 
Sci 17(6) pp 1169-1220 (2006). 



[14] A. Ambainis, Quantum walk algorithm for element distinctness, SIAM 
Journal on Computing, 37:210-239, 2007. Also arXiv:quant-ph/0311001. 

[15] A. M. Childs, Universal computation by quantum walk, Phys. Rev. Lett. 
102, 180501 (2009). 

[16] S. E. Venegas-Andraca, Quantum walks: a comprehensive review, Quan¬ 
tum Information Processing vol. 11(5), pp. 1015-1106 (2012). 


